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ABSTRACT 


The stability of plane Poiseuille flow was studied 
using theory developed by Harrison. A similarity trans- 
formation was introduced which reduces computation time 
and provides better insight into the basic relations. 
The stability of the flow was examined from a Lagrangian 
viewpoint. Instability was found to be progressive in 
nature and three distinct levels were identified, namely 
incipient, critical, and fully developed instability. 

Results show that the critical Reynolds number can be 
lowered indefinitely if certain types of perturbations 
occur. Specifically these involve relatively abrupt 
changes in amplitude. This provides a possible expla- 
nation for the disagreement between earlier theory and 


experiment. 
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All quantities are expressed in dimensionless form by 


the use of a natural system of consistent units in which 


channel semi-height is the unit of length, the volumetric 


mean velocity of the fluid is the unit of velocity, and 


the density of the fluid is the unit of density. Then 


all other consistent derived units are fixed accordingly. 


A* 
e 


G 


wave number amplitude defined in Eq. 2.4 
2.71828...base of the natural logarithms. 
Stability parameter defined in Eq. 1-30. 
complex stream functions. 

transformed stream functions. 

complex auxiliary function defined in Eq. 2-14. 


att? 


unit vectors along x, y, and z axes, 
respectively. 


the imaginary unit. 


complex auxiliary function defined in Eq. 2-15. 


transformed auxiliary function defined in 
Gree 75 


Reynolds number based on volumetric mean 
velocity and channel semi-height. 


transformed Reynolds number defined in 
Gy we— 8s, 


complex auxiliary function defined in 
Bae. 2- | Se 


time, 


flow velocity. 
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COMplcreVeCuOrepOcential OF perturbation 
flow defined in Eq. 3-2. 


coordinates in direction of mean flow, 
normal to walls and transverse to the mean 
flow, respectively. 


coordinates in moving reference frame. 


complex wave number of the perturbation in 
Cede neerron . 


transformed wave number defined in Eq. 2-3. 


complex wave number of the perturbation in 
the z direction. 


complex frequency of the perturbation in 
a fixed reference frame. 


complex frequency of the perturbation in the 
moving reference frame. 


transformed complex frequency defined in 
Eee 2d, 


phase angle parameter defined in Eq. 2-12. 
amplitude parameter defined in Eq. 2-19, 


angle of plane of perturbation with 
Lespect CO Ny plane. 


angle of resultant growth wave number vector 


he With respect to X axis. 


angle of oscillation wave number vector XY 
with respect to x axis. 


growth wave number vector. 

oscillation wave number vector. 

wave number phase angle defined in Eq. 2-9. 
wave number phase angle defined in Eq. 2-7. 


wave number phase angle defined in Eq. 2-2. 





I. THEORETICAL BACKGROUND AND APPROACH 


A. BACKGROUND 

This research deals with the instability of plane 
Poiseuille flow, that is, plane flow between infinite 
parallel plates. The mean velocity of this flow is given 


by the expression 
3 2 
U = x(1-y ) : (1-1) 


The stability of such a flow field is determined by 
Superimposing upon it an appropriate perturbation and 
determining whether this perturbation tends to grow or 
decay over time. In the present case the perturbations 
are expressed by a complex vector potential which is taken 


to be of the form 
W = [JG(y) + KH(y)] exp(ax + 8z + yt) . ce 


The complex constants a and 8 fix the spatial charac- 
teristics of the perturbation and may be arbitrarily pre- 
scribed whereas the complex constant y fixes the response in 
time and must be found by solving the vorticity transport 
equation. Moreover, since a, 8 and y are all complex, they 
can be resolved into real and imaginary components in the 


form 


10 


ne oe 








ee 


R I 
8 = Bp + iB, (1-4) 
eee lYy (1-5) 


Thus the spatial characteristics of the perturbation 
are seen to be completely defined by the four constants 
Ap, Ar; BR; By . In addition, the mean flow is charac- 
terized by its Reynolds number Re. 

Harrison's original analysis [Ref. 1] showed that the 
perturbation growth rate in time as seen by a fixed observer, 
and as expressed by the parameter YR» is a definite function 
of the five parameters Op, A, Bp, By and Re, which charac- 


terize the perturbations and the flow. Thus 


oe - Yelp s Te Bs By» Re ] ° (i) 


femee ROGRESSIVE INSTABILITY 

In a further development of Harrison's original approach, 
Section II of this thesis shows that three significant levels 
of instability can be defined which are termed incipient, 
critical and fully developed instability. The definition of 
these terms depends, in part, on the algebraic sign of Op: 
However, Harrison showed that negative values of Op have a 


definitely destabilizing effect. Consequently the present 


om 





analysis is restricted to the critical case of negative Ap 
For this case the three levels of instability correspond to 


the following levels of Yr: namely 
Incipient Instability Yq = 0 (1-7) 
Cratiacal Instability YaJe = -Op (1-8) 


Healy ewe oped Eee. 
Instability YRJp oe a 


Numerical solution of the vorticity transport equations 
enables us to find the three corresponding Reynolds numbers 


at which the above stability levels are reached. Thus 


Incipient Instability Re); = RE; (ap, a7, Bp, By] (1-10) 
Getetecal Instability Re) ¢ = Rec lap, a7, Bp, By] (1-1) 
Fully Developed _ - 

Instability SS eee Ree rote) 


C. TRANSFORMATION OF PARAMETERS 
miecection Ll of this thesis, a transformation is 


developed which relates the original parameters Ap, Ars Bp» 
Br» Re and Yp to an alternative set of parameters which are 


* 


* * * 
Symbolized as A , # , 8, Re and Yp: This transformation 


eZ 





can be expressed in several alternative but equivalent forms. 


In the present context it is convenient to write 


* * 
Op = A cos(f +8) (T=13) 
* * 
a7 = A sin(g +6) (1-14) 
2_aA2 a) sea z 2 ‘ 
Bp = a fel (Gleaher te 4nG = S Teo tcosiZ0-kK COS2())) to) je (raseoy 
2_aA2 2 2. Die * 2 i 
By = ee «(tL (1 ie tees ol -cos2f +k cos2(p +6)} (1-16) 
Re = 2s ee (ie) 
and 
k 
ee KY - (1-18) 


The important fact about this new set of parameters, 
which are somewhat loosely termed the starred parameters, 1s 
that their use permits the fundamental vorticity transport 
equation to be simplified. Specifically, the relation 


analogous to Eq. 1-1 reduces to 


k kok * x 
Yp = YplA ,ae , 8, Re |] (1-19) 


The relations analogous to Eqs. 1-7, 1-8, and 1-9 reduce 


EO 


les 





* 
Incipient Instability YpJq = 0 


; * Z * 
ere ical Instability Yale [ae 
nuely Developed 2 —. 
Instability YR’D 2 
Finally, the relations analogous to 


Simplify to 


* 
Incipient Instability Re, = 
4: ee * 
Crnutcatelnstaurlaty Reo = 
Fuliyeveve loped nist 2 
ims tab 1 i aty D 


The remarkable feature of Eqs. 


(1-20) 


(ieee <a (ail) 


bol 


3% x 
= z A Coste Feuer) 


Eas. 1-5, 1-6, and 1-7 


RE;[A , § , 8] Gs) 
* * * 
= Reo [A ’ p ’ 6 | (diaz) 
* * * 
Ren [A ’ p ’ 6 | (1-25) 


Loe taeOugn eh 2 oes 


that none of these relations involve the parameter xk. 


Thus the number of independent parameters has been reduced 


PeemecOur in Eqs. 1-10, 1-11, and 1- 


1-24, and 1-25. 


12%te three in eqs. 1-22. 


This represents a very significant simpli- 


fication of the problem, especially in view of the tremendous 


computational burden which these equations involve. 
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D. PHYSICAL SIGNIFICANCE OF STABILITY BOUNDARIES 

The stability boundaries Re), Re )a, and Re) 
symbolized by Eqs. 1-23, 1-24, and 1-25 have physical 
Significance which can be interpreted in a straightforward 
manner. Eq. 1-17 shows that Re. can be regarded as the value 
of Re which corresponds to the reference case kK = 1. Thus 
Re ) 7, for example, is the Reynolds number of incipient 
instability for a perturbation which is characterized by 
the given values of parameters A. O°. and 6 and by the 
reference value K = 1. Since the transformed quantities 
A, a. and 6 may, at first, seem to be somewhat abstract 
in character, it is helpful to go back to Eqs. 1-13, 1-14, 
1-15, and 1-16 and ascertain the corresponding values of 
the original untransformed parameters Ap, G7, ne By - 

However, there is far more to this solution than just 
the above reference case, K = 1. In this connection, 
Meme i-i7, when taken in conjunction with Eqs. 1-23, 1-24, 
and 1-25, reveals a most important result. It shows that 


* * 
for given values of parameters A , @ and 8, and hence 


for the corresponding values of Re )y, Re), or Re), 
the corresponding actual Reynolds numbers Re)y; Rea, or 
Re) p can be lowered indefinitely, simply by increasing 
parameter kK to any desired extent. Notice that such a 


shift of the stability boundaries, while it involves no 


* * 
change in parameters A , § , and 8, does involve changes 


ILS 





in Ap, A> Bp: and By. It is therefore important to 
Summarize the nature of these changes in the clearest 
possible way. 

For this purpose, it is useful to regroup the four 
quantities ap, %,, Bp, and 8; into two vectors Xp and Ay 


defined as follows: 


AR = Lop + KBp (1-26) 
xy = Lo, i KB, (au27) 


fcieand k are unit vectors in the x and z directions, 
respectively.) 

Clearly Ap represents the Spatial growth rate in vector 
form, that is, in terms of magnitude and direction, while 
AY represents the spatial oscillation rate in like terms. 
Each of these vectors is characterized by a magnitude and 


a direction. In this case the magnitudes AX, and dy turn 


R 
out to be governed by the relations 








eZ 
Ag = “ {(((1-K2) 2+4«2%sin29)1/2- (1-«2))+2c0876'} (1-28) 
2_ Av? DO? ene 2 2.8 
Ay ae, Gl Ol yendc=sain 6 | -(l1-K"))+2sin’# } . (1-29) 


R and Ay Which the 


above vectors make with respect to the x axis turn out to 


Likewise the two corresponding angles A 
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— a 


oe ae a, 


a re 


7 a > 
oo ‘_ eae ; - - 7 a 





_ - 
, oe 
: a 
; 
a 
7 
- . 
7 7 - : 
- 7 oe 5 - - 
- : ew 7 od - 7 / 
Dy ie Gs | 
~ a 
| ™9 - - 7 
7 i 
| SS 4, ts 
7 . 2 7 — 
| . <a 
| > i = 
| 7 7 - ~- ; 
OAae - —s . i 7 - 
- 7 7 
7 7 C y 7 © 7 > - 
7 
7 7 / : 7 ae ; 7 ; 
7 7 > a : / 
- - 7 - an 7 - 7 Oo ~~ 7 oe a oe a 
a 7 7 7 YT as ~~ RS - - ee a a oe 
a 7 —~ mee _ on Beene se 86h else CUD te; 


be governed by the relations 


Oe Die eae leer e 2 * 
tan“A, meiehek ) t4eesin 0] MeEGo2 be -k COSZ(f +) (1-30) 


eel + eos =o) | 


es 2 weg 2 ely ee 2 - 
ee PeGlaiceeudi San 7 0) -cos2$ +k cos2(g +86) (1-31) 
K {1 - cos2(@ +86) ] 
Thus the spatial form of the perturbations is now fully 
characterized by the four transformed parameters ps hyo Ap; 


and Ay which are in some respects more convenient than the 


four original parameters Ap, 7, BR; and Br. 


E. EFFECT OF VARYING PARAMETERS 
The connection between the above perturbation charac- 
teristics and the Reynolds number is still expressed by the 


relation 


+ 





Re (1-32) 


Ke 
It is very instructive to study the trends revealed by 
Eqs. 1-28 through 1-32 when parameters Re a. and @ are 
held constant while kK is allowed to increase. Eq. 1-32 
reveals that Re)y; Re)o, and Re) p can be decreased indefi- 
nitely in this manner. On the other hand, Eq. 1-28 reveals 
that any such decrease in Reynolds number always entails a 


corresponding increase in the quantity A Recall that A 


R° R 


ile? 





represents exponential growth rate in space. This says 
then that stability boundaries are not absolute in character 
but depend significantly on the "abruptness" of the pertur- 


bation in space as measured by parameter i The greater 


R° 
this abruptness parameter, the lower the Reynolds number 
at which instability can occur. 

Conversely, if the permissible magnitude of Ap be 
limited in some definite manner, the reduction that can be 


achieved in Re)y> Re) and Re) 5 Will be correspondingly 


G 
limited as well. In that case, a systematic exploration 
over appropriate ranges of the parameters Ae O° and 6 
Should ultimately reveal corresponding ultimate stability 
limits and, in particular, some critical Reynolds number 
below which ~ instabilities occur. Notice, however, that 
such a critical Reynolds number is never absolute, but is 
always contingent upon the restriction that has been placed 
upon parameter ps 
1. Restrictions on hp 

The most obvious and direct restriction that can 
be placed on hp 1s Simply to limit it to some fixed value 
or, for study and comparison purposes, to some series of 
Successive fixed values. In general, the boundaries which 
correspond to incipient, critical and fully developed 
instability will then depend on the designated value of 
Ap: The higher this value, the lower the values of Re 


at which the above boundaries will occur. 
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Any such restriction of the magnitude of A, implies a 


R 
corresponding restriction on kK. To show this, invert 
Eq. 1-28, solving for k as a function of Ap. The result 


1s 


- psy’ - cos*# ps)’ + sin?g" 


(=) - cos“f 5 sin?) 
A 


K (pliq Soy) 


This relation may be used in connection with Eq. 1-32 
to express the final Reynolds number at which the designated 
Stability boundary is reached. For incipient instability, 


for example, this boundary may be expressed in the form 


* * * 
Re; (A ,f 8) (1-34) 





Analogous expressions apply to the boundaries for critical 
and fully developed instability. 

Equation 1-34 shows quite clearly that the stability 
condition in question depends on the three characteristic 
parameters A, A, 8 of the perturbation as well as on the 


limiting value assigned to parameter hp: This procedure 


LS 











of calculating stability boundaries for various assumed 
combinations of ao A. 8 and Ap has been carried out for 
several typical cases and the detailed results are sum- 
marized in Section IV of this paper. Of course, these 
examples, while representative, merely scratch the surface 
Samene stability problem. The complication remains that 
the true and ultimate stability boundary represents the 
lowest possible Reynolds number at which an instability 
can just occur. This implies that all possible combinations 
of parameters ie g° and 8 must be examined to determine 
the particular combination which, for a given limit on p> 
yields a stability boundary at the lowest possible value 
of Re. In other words, the true stability boundary amounts 
to the envelope of all the individual stability boundaries. 
Each individual boundary is characterized by some 
Specified combination of A’, a” and § and, of course, also 
Ort pe: Since there is an unlimited number of such combi- 
nations, the amount of calculation involved in establishing 
the desired stability envelope is prodigious. Needless to 
say, no such attempt was made in the present thesis to 


accomplish anything so ambitious. 


PoC OPE OF PRESENT RESEARCH 

The present research was restricted to the more modest 
and realistic aim of calculating stability boundaries for 
a few specific and typical combinations of An o. 8 and 


ps This goal has been successfully attained. 
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G. PARAMETER G 
1. Definition of Parameter 

A detailed study of the relations summarized by 
Eqs. 1-28 through 1-32 reveals the possibility of express- 
ing a restriction on the permissible magnitude of p in a 
rather subtle and indirect way, through a change of variable. 
The particular algebraic form which the above relations 
assume suggests the utility of defining a new parameter, 


called G, as follows: 


@? = 1/2 {[(1-«2)* + 4n*sin°e}1/4 - (1-«*)} Caan 


This relation can be readily inverted to give 


Mame (Goes 
Ge Sim 8 


fe Utilization of Parameter G 
Equations 1-35 and 1-36 may be used to eliminate 
parameter « from Eqs. 1-28 through 1-32, replacing it by 
the new parameter G. In this way the following results are 
obtained. 
The vector amplitudes Ap and dy turn out to be 


related to the new parameter G in a fairly simple fashion. 


Zo 





: 
; 


a 





Thus 


> 
il 


A ((G2 + cos“g }]2/2 (eam 


hace + ee ya (ee Sie) 


> 
tl 


On the other hand, the angles An and Ay 


Simplified by the use of parameter G. Fortunately these 


ayesnot 


quantities are less significant than the preceding ones. 


The governing equations become 


i Die meee 2. ee eae 
tan“A, BliGenoineo(G -Sima) JaGa(Gs ti) sim: (eco ) (1-39) 
Geta t1jeos =p eu) 
and 
ae. 2 2 Lee ce. 2 Ce 
tan“A, (Gt sineoGe cos 0) ) tGe (Get iicos aGl. 70) (1-40) 


G (Getlisin ( +0) 


The important Reynolds number relation below is once 
again simple. For definiteness it is written specifically 
for the case of incipient instability, by analogy with 
Eq. 1-34. Similar expressions apply also to the boundaries 


of critical and fully developed instability. Thus 


Ze 


Re), = G"_+ sin’e Re, (Af ,8) (1-41) 
Cease: 

Equation 1-41 shows that the stability depends on 
the particular parameters io o. 8 and on the limiting 
value assigned to G. Hence G plays a similar role in 
relation to Eq. 1-41 that hp plays in relation to Eq. 1-34. 

The results summarized elsewhere in this thesis 
are presented primarily from the perspective expressed by 
Eq. 1-34. The alternative version shown by Eq. 1-41 is 
included in this discussion because of its theoretical 
interest, but this version 1S not used in the presentation 
of calculated results. 

Notice that in either version, assuming some assigned 
limit for Ap or G, an exploration is still required over the 
domain of parameters AY, g” aniags tOuLind the pauticulan 
combination that yields the minimum Reynolds number. Of 
course, such extensive exploration could not be undertaken 


im the present thesis owing to time limitations. 


eee e DUCTION TO CLASSICAL THEORY 
It 1S pertinent to note that the classical theory of 
the stability of plane Poiseuille flow amounts to a special 


case of the more general theory discussed above. It 
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emeunts, in fact, to the special case for which 
eee. (1-42) 


Study of Eq. 1-28 reveals that Eq. 1-42 can be satisfied 


@emand Only if we set 


e= 0. i=43) 


and 
Pozo. (1-44) 
From Eqs. 1-28, 1-32, and 1-42 we may infer also that 


ee | Ct (1-45) 


Moreover, we also find under these conditions that 


eee 0 4 and@a, = 0). (1-46) 


It is evident that the general theory discussed in this 
thesis is immensely more comprehensive than the classical 


theory as limited by Eqs. 1-42 through 1-46. 
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II, SIMILARITY TRANSFORMATION 


The perturbation characteristics are fully defined, as 
in Ref. 1, by the four real parameters Ap, Ay, Bp» and By: 
Ap and a, are the components of the complex wave number of 
the perturbation in the x direetion and BR and By are the 
components of the complex wave number of the perturbation 
mimeene Z direction. 

The above parameters satisfy the following relations: 


 . 
a Ser 2 Segawa (2235) 
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~ 
@ 
I 


: ; * * 
A very useful alternative set of parameters is Ap» Ars Q 


aad (°K . 


eZ 
a 1s defined by 


G, = 0? + a? : (2-3) 
* * * * 
am. oO , Op and a; are defined by 
* & A” : - 
a 
a =Ae = Ap + lay (2-4) 
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Op and a, are the components of a , the transformed complex 
wave number parameter and kK is the transformed perturbation 
amplitude parameter. 


Kk and @ are defined by 
oF cca. (2) 


An and K are positive by definition. 

Given the original parameters Gps Ar; Bp and By: the 
transformed parameters an ays 8 and k may be deduced from 
equations 2-1 through 2-5 as follows: 


2 1/4 
Z Z Z Zi 
Av = [(ap - af + Be - 84) + 4(apa, + 828,)7] (2-6) 


Zio + Bae) 
o° = vl/ 2 sareean —— (2-7) 


KAW = [ag + atyll? (2-8) 
1 
iy = srevan( | (259) 
* * * 
Ap = A cos @ (2-10) 
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Q 
ll 


te x 
A sin @ (iz gis) 


(p - 6°) (2-12) 


cD 
Il 


The governing equations of Ref. 1, using the five 
independent parameters, Re, Op,» A, Bp» By are as follows: 


Using the auxiliary expressions below, 


2 2 
THe Se .- apti-y*) (2-13) 
at + gt 
Oo = Vee Ley) (2-14) 
J(y) = (a + 87)T(y) - 3a (2-15) 


the fundamental vorticity equation becomes 
1 yiv tt '" Z Z z 
eee * ty) RY + S(y)B] - y[R" + (2° + BDH] = 0. (2-16) 
The associated vorticity equation is 
d. ? = B i tt ' 
ie" «| (UCT (y)-y)G = wrage Re HN! +(T(y)-y)H'-3ayH] . (2-17) 
GQ + 


The number of independent parameters in equations 2-13 


through 2-17 can be reduced to four by utilizing the 


oy 





following functional transformations. 


H(y) = « te7?°H* cy) 


= 


Gigjueun oe + Eero 


plier cucion of Eqs. 2-18 through 2-23 
Sommeron, Eqs. 
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«Lobes 
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* 
Daye) 


O 
= ae 


Re 


* * 
Ly) eye) 
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x" (y) 
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(1-y*) 


a * * 4 
a T (y) - 5a ehh 


(2-18) 


(2.2 10S) 


(2-40)) 


(2-21) 


(Za) 


(2-23) 


imto the emerat 


2-13 through 2-17 yields the three auxiliary 


(2-24) 


(2-775) 


(2-26) 
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The principal vorticity transport equation becomes 


V ; 2 


l real * #11 a 3) *& *! * * 
fae CU tC U(yJH 6+ Ch (y)H J] - y (H ta HJ = 0 C2527) 
Re 


The associated vorticity transport equation becomes 


l at * * 1 1 «i <e 
[—-G +(T (y)-Yy)G = — eH logy cl 
Re A Re 


* *& 
seen yi ae (2-28 


Equations 2-24 through 2-28 now involve only four 
* * * 
parameters Re , Gps Oy and 6, «, the fifth parameter, has 
Gameelled out. Kk becomes part of the solution again during 


the reverse transformation of results from starred parameters 


to the original parameters. 


Zo 





Eli. oAerL PieeeeR TT ERION 


A. BACKGROUND 

Studies of the stability of Poiseuille flow have used 
various criteria for determining the stability of the flow 
from the solutions obtained. The growth rate in time, Yp> 
is usually used when there is no real-exponential spatial 
variation [Salven and Grosch, 1972]. When exponential 
growth in space has been included [Garg and Rouleau, 1971] 
the real part of the spatial wave number has been used to 
give the instability but this procedure, while seemingly 
plausible at first inspection, cannot be really justified 
with any rigor. The Lagrangian approach described below 
is believed to be a Superior method for dealing with this 
case. In other cases, stability has arbitrarily been 
evaluated with respect to a frame of reference moving down- 
Stream at the phase velocity of the perturbation but again, 
this procedure has no strict rational justification, and 
especially so in connection with the fully three-dimensional 


perturbations considered in this thesis. 


B. LAGRANGIAN REFERENCE 

For perturbations that are both oscillatory and have 
exponential rates of growth or decay in the streamwise and 
transverse directions a Lagrangian frame of reference 


proves useful. The fluid particles have velocities varying 
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meee to «1.5 depending on their distance, y, from the 


walls. The velocity distreunueron iS given by 


U = (1-y*) ean 


Consider a coordinate system moving in the x direction 
With the mean velocity of a given fluid particle. Let y 
be the mean vertical coordinate of the moving particle. 
Then the velocity of the moving reference frame is the same 
as the velocity of the streamline along which the above 
Beaterrcle moves and 1s given by Eq. 3-1. Let x', y, z, t 
Dewtne coordinates and a, 8, y' the complex wave numbers with 
respect to the moxing axes. The form of the perturbation 
vector potential for a given eigenvalue obtained as a 


solution is 

W = jG(y) 1p kH(y) expe + Bz + Vt). , (3-2) 
iiemweomplex frequency y' seen from this moving reference 
is different than from a fixed frame. To relate y' to y the 


perturbation vector potential is written in the moving 


frame and transformed into the form for the fixed frame. 
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W = (jG + kH) exp (ax! sey ete tee) 
= (jG + kH)exp(a(x-Ut) + 82 + y't) 
= (jG + kH)exp(ax + 8z + (y' - aU)t) 
= (jG + kH) exp (ax ery ce) (550) 
Therefore 
1 eee (3-4) 


solving for y' and splitting into real and imaginary 


Pamts yields 


tan (5-5) 


Yt = ¥; 7 a,U (3-6) 


if YR 1S positive, zero, or negative, the perturbation is 

Said to be unstable, neutral, or stable, respectively, with 
respect to the moving reference frame. Thus, the value of 
YR is taken to be a measure of the stability of each eigen- 


value obtained. 
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fee eRANSFORMED STABILITY CRITERION 


Consider now the transformation to the starred parameters. 


The condition of stability is determined by the value of YR 


which Eq. 3-5 gives as 


! = 
oo pee pe 
Now 

t 

7. 

eee 

_ * 

Yeo SR 

a * 

a = ae 


Peberrtution of Eqs. 3-7 into Eq. 3-5 yields 


* 
Op is defined by Eq. 2-10 as 


k & < 
dp = A GOs (0 on 


Therefore 


ni 
YR 


* * * 
ie A cos(g +6) 


(3-5) 


(3-7) 


(3-8) 


(2-10) 


(3-9) 





ie three stability boundaries, incipient, critical, 
and fully developed are determined by the condition that 
! 


@xmists when i 0. Setting Eq. 3-9 equal to zero and 


* 
solving for YR yields 


& & sk 
i UA cosi@pe +9) 4 (3-10) 


Now incipient instability corresponds to zero growth rate 
With respect to a coordinate system which moves with the 


flow velocity at the wall, which is zero. 
* 
YpJz = 9 Cap 


Critical instability corresponds to zero growth rate with 
respect to a coordinate system which moves with the mean 


Wemocity of the flow, which is unity. 
- ~ ~ 
YRI¢ = -A cos(f +8) (3-12) 
Fully developed instability corresponds to zero growth 
rate with respect to a coordinate system which moves with 
the velocity of the flow on the centerline. 
* 3 *% * 
Yap = - zA cos(f +86) (sS)3)) 


2 


34 


ie kesules 


A. TRANSFORMED PARAMETERS 

Equations 2-24 through 2-28 were solved on the IBM-360 
computer and the most unstable growth rate, Ye for a given 
ye, 0, Re, AY was Gbtained, ihe botmdaries for inciprent 
and critical instability were determined by the criteria 
cupbaimed in Section [I. Fully developed instability did 
not occur for the case studies. Two values of a" were 


a “ 
Studied at various values of 9, Re and A, 


ie 90° 

Values of 8 explored were 0, 1, 2, 3, 4, 5, and 6°. 
Because one-degree increments of 8 yield graphical results 
that are extremely cluttered, this study will present the 
results only for 6 = 0, 3, and 6°. Figure 4-1 shows the 
Stability boundaries for a” = 90°. Thevetfect of changemye 
8, While holding a” constant, can be seen. An increase in 
8 causes a degrease in ee for both ineipient and critical 
instability. Also, for a given 9, the boundary for critical 
instability occurs at a higher Re than does the boundary 
for incipient instability. 

There is only one curve for 6 = 0°. In this case 
the criteria for incipient, critical, and fully developed 


instability turn out to be identical. Equation 3-10 shows 


* * * 
oe eee cosip 0) (3-10) 
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emeecos(90 + 0°) = 0; thus the criterion for stability on 
the three boundaries 1S Yr = 0. 

Note that the boundary of incipient instability for 
86 = 3° shows an abrupt distoncinuity represented by segment 
ab. This discontinuity iS similar to that obtained by 
Harrison. It is expected that extension of the incipient 
boundary for other values of 8 would reveal the same 
characteristic. 

Figures 4-2 and 4-3 are plots of the growth rate, 
Vp? versus Re for 6 = 3° and 6°, respectively. Both plots 
show the locus of points that represent the boundaries of 
nena ent ao €YGitical instability. Let can be Seemmenan 
fully developed instability is not reached even at 
Re = 100,000. 

2. pf = 95° 

Due to a lack of time only two values of 6 were 
explored, 8 = 0 and 3°. A comparison of Fig. 4-4 with 4-1 
Shows that the stability contours follow much the same 
pattern for both cases. However, increasing g” to 95° 
causes a corresponding decrease in Re 

Figures 4-5 and 4-6 show the growth rate as a 
Pumetion of Re for 8 = 0 and 3°, respectively. The locus 
of points that represent the boundaries of incipient and 
critical instability can again be seen. Fully developed 


* 
instability is not reached at Re = 100,000. 
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B. RESULTS TRANSFORMED TO UNSTARRED PARAMETERS 

In order to present the results in a manner consistent 
with Ref. 1, it is necessary to transform the results from 
starred to unstarred parameters. This transformation puts 
the results in a more easily understandable form. 


From Fig. 4-7 the following relations can be deduced. 


r 
(—*) = 1/2(K* + [(1-«7)* + 4x?sin?e] 747+ cos2p") (4-1) 
A 


2 
A 
(2) = W/ 2c r ialee + ieosinee ye = cos2f) (4-2) 
A 
tan“Ap 7 Clee sence teas Wee cos7(t Fe) (es 


k(t > cosZy ed 


2 2 


Dix 2 Wye 2 ee - 
tan?A, Sa Love fds sine s) | / -cos2@ +K cos2(@ +6) (4-4) 


K (l= eosZ Glee) ] 
Ap and Ap are the magnitude and direction, respectively, 


of the perturbation and growth vector. A, and Ay are the 


I 
magnitude and direction, respectively, of the perturbation 
wscillation vector. 

1. Perturbation Rate Vectors, Magnitude 


* 
Figure 4-8 shows G as a function of Re/Re with 8 


aS an independent parameter; the curves are valid for all 
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Figure 4.7. Vector Diagram for Parameter Transformation 
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values of a. For any fixed value of 8, G decreases as 
the Reynolds number ratio increases. In the case of 9 = 0° 
G decreases with increasing Reynolds number ratio until 
pene” = 1. G then maintains a constant value of zero for 
further increases in the Reynolds number ratio. Although 
Values of @ above 6° were not explored, values of 8 up to 
60° are shown here to demonstrate the trend as 6 increases. 
2. Perturbation Growth Rate Vectors, Direction 
Pie-quantitaes of tan“Ap and tan“A, are fixed by 
Eqs. 4-2 and 4-3, respectively. However, if tan“A, be 
Specified, this does not fix Ap uniquely as there are 
four angles, one in each quadrant, which have the specified 
value of tan“A,. Similar considerations apply also to the 
other angle Ay. Hence, to determine Ap and Ay uniquely 
it is necessary to consult auxiliary relations which fix 


the quadrant in which these angles really fall. 


The components of Ap Which fix angle Ap are 
* * 
a, = ApcosA, = KA cos(f +6) (4-5) 
and 
: * 
Bp : ApSinh, = GA cosy . (4-6) 


The components of At which fix angle Ay are 


* * 
a; = A cosh, = KA sin(@ +6) (4-7) 
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and 


° * ° 
B, = A;SinA, = OA sin . (4-8) 


* 
Moreover, in this study, the angle (§ +8) has been 


restricted to lie in the second quadrant so that 


{ A 
Oo 


eel) 


Iv 
oO 


ay (4-10) 
Recall that negative values of ad, were shown by 
Harrison to be destabilizing. That is why the present 
Study is restricted to negative values of Op. 
In order to determine the algebraic signs of 
components Br and Brs 1t is necessary to bracket the range 
of the angle yw. A study of Fig. 4-7 reveals that for 


positive values of 8, the following limits apply. 


a, (ae 
ra 0 
Cone ; (4-12) 


K +» 0 


It is also evident that the x-y plane is a plane of 


Symmetry and that therefore a reversal of the perturbation 
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characteristics with respect to the z axis is permissible 
and leaves the essential features of the solution other- 
wise unchanged. This amounts to saying that the angle wy 
can always be changed by +180°, with no significant effect 
upon the solution except for a reversal of the perturbations 
with respect to the plane of symmetry. For definiteness 
in this discussion, however, we limit the angle w as inci - 
cated by Eqs. 4-9 and 4-10. It is then evident that aite 
above auxiliary relations, along with the basic relations 
Geemegs. 4-3 and 4-@ , now suffice to f£1x An and Ay uniquely 
for any assigned values of the parameters A”. 8 and kK. 
Figure 4-9 shows A, as a function of Rene EOF 
a = 90° with 8 as an independent parameter. This plot 
shows a constant value of A, = 90°, fom 9 = 0° and Re < Re. 
When Re > Re. the vector magnitude, Ap 1s"’zgemwo. Gor all 
Other values of 8 the perturbation growth rate vector, dev 
rotates from near the transverse to near the upstream direc- 
tion as payne ™ increases. 
Figure 4-10 shows the rotation of the perturbation 
growth rate vector as Re/Re- increases, for g” = 95°, Note 


* 
that when § = 90°, A, varies between 90° and 180° whereas 


R 
* 
when § = 90°, A, varies between 90° and 270°. 
5. Perturbation Oscillation Rate Vectors, Direction 
Figures 4-11 and 4-12 are similar plots showing 


the rotation of the oscillation rate vector with changing 


* * 
Re/Re for §@ = 90° and 95°, respectively. Comparing 
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Fig. 4-11 with 4-9, for f° = 90° and ® = 0° both A, and A, 
have constant values when Devine < unity. 
4. Stability Boundaries in Unstarred Parameters 

Figures 4-13 and 4-14 show the stability boundaries 
for unstarred parameters. Ap HSeee05 for bern plots. 4 
comparison with Figs. 4-1 and 4-4 shows that the character 
of the boundaries remains unchanged. However, Reynolds 
number is greater than starred Reynolds number. 

Although it is not shown here, it was found that 
increasing Ap causes the boundaries to move to the left 


and up. In other words, an increase in Ap causes an 


increase in dy and a decrease in Reynolds number. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The research reported here is based largely on the 
theory developed by Harrison in Ref. 1. However, this 
theory is further extended in the present work by the 
introduction of a useful similarity transformation. The 
transformation reduces the number of independent parameters 
required in the computer solution from five to four and 
thereby significantly reduces computation time. 

In Ref. 1, Harrison showed that negative values of Op 
are destabilizing. The results presented here show that, 
for negative values of Ops the critical Reynolds number for 
plane Poiseuille flow can be lowered indefinitely by 
proper selection of perturbation characteristics, provided 
that the corresponding increase in p 1s acceptable. Thus 
the stability boundaries are not absolute in character but 
depend significantly on the "abruptness" of the perturbation 
in space as measured by parameter Ap: The lowering of the 
critical Reynolds number in this way provides one possible 
explanation for the disagreement between earlier theory 
and experiment. 

The stability of the flow along a particular streamline 
has been shown to depend on its velocity. Negative values 
Or: Op yield the greater instabilities and streamlines with 
the lowest velocities, those nearest the walls, will be 


the most unstable. This seems to agree with experiment. 
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Peeceonaing to Sehlichting, transition from laminar to 
turbulent flow is "characterized by an amplification of 
the initial disturbances and by the appearance of self- 
Sustaining flashes which erpnieine from fluid layers near 
the wall along the tube." 

Instability was found to be progressive in nature and 
two of the three defined stability boundaries were located, 
incipient and critical. Further research needs to be done 
for a wider range of parameters A f°, and @ to find the 
combinations that correspond to minimum Reynolds numbers at 


the above stability levels. 
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APPENDIX A 


USE OF THE COMPUTER PROGRAM 


It was found extremely useful to precompile the program 
on a disk thereby avoiding the inconveniences of reading in 
the complete card deck for each run. An additional advantage 
1S a reduction in turn-around time of considerable magnitude. 
The following will give procedures and hints that can be 
found in the W. R. Church Computer Center but require time- 
consuming search. 

1. Pre-compiling Program 

To compile the program, the following was read 
into the system 
// Green Job Card,Time =(0,59) 
eee xcC FORTCL 
poe ORT.SYSIN DD * 
/* 

Program Card Deck goes here. (no data) 
Pe aeNK.. SYSUMOD DD DSNAME-S2593.LIB(POIS) ,DISP=(NEW, KEEP) 
// UNIT=2321,VOLUME=SER=CEL006 , LABEL=RETPD=220, 
ee OP ACE=(CYL,(6,1,1) ,RLSE) < 
/* 
2.  Erovran Execution 

Once the program is compiled, running the program 

consists of punching data cards in the namelist format and 


reading them in with the following deck of cards. 
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// Green Job Card,Time=(0,59) 

// GO EXEC PGM=POIS,REGION=178K 

//STEPLIB DD DSNAME=S2593.LIB,DISP=SHR, 

// VOLUME=SER=CEL006,UNIT=2321 

//FTO6FO01 DD sysout=A, DCB=(RECFM=FBA,LRECL=133,blksize=3325S) 
//FTOSFOO1 DD * 

LIST N=30,REY=3000, TH=.052360,ASTAR=1.8,PHIS=95, SEND 

/* 


Note: Column 1 is blank on the list card. 


Three decks of these cards were used with each having a 
different job name, i.e., NEWBY 64A, NEWBY 64B, and NEWBY 64C. 
This proved useful as three jobs could be loaded at one time 
and one could keep track of what was already printed and 
what remained to be processed. It was also found to be 
useful to have three sets of job cards with each set having 
a different time. For Time=(0,59), S59 sec., one list card 
(data) was inserted. This was used for quick turn-around 
time and only a few points were being explored. For 
Time=(2,00), 2 minutes, three list cards could be read in 
and for Time=(4,00) six list cards could be used. Occasion- 
ally the four-minute time parameter would terminate execu- 
tion after five list cards had been processed. 

3. Program Alteration after Compilation 

If changes were to be made in the program the file 


was scratched and a new file established with the changes 
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Miconporated. To scratch the program on file the following 


deck was used. 


// Green Job Card 
// EXEC PGM=IEHPROGM 
//SYSPRINT DD SYSOUT=A 
//DD1 DD UNIT=2321,VOL=SER=CEL006,DISP=SHR 
ZoroiN DD * 
SCRATCH VOL=2321=CEL006 , DSNAME=S2593.LIB,PURGE 
/* 
Note: The scratch card begins in column 35. 


In all three decks the name of the program (POIS) 
appears. The choice of a program name is an individual 
choice but once chosen it must appear the same in all 
card decks. The only other item that appears with unique- 
ness is the individual user number. In the context of this 
paper that number was 2593 and must agree with the user 
number on the job card. 

There are two possible selections on input parameters 
for obtaining data. Both are used in this study. It is 
Ret and vary ie to 


* 
possible to select values of § Q 


3 3 3 


find a point. This method was used first and worked well 
when obtaining a solution from the computer. The problem 
arises when interpreting and transforming the results. It 
proves useful to have values for fixed he and this method 


dees MOt provide this easily. 
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* 


* 
An alternative technique is to select 9 6, A Weanid 


3 3 


* 
then vary Re . To construct Figures 3-1 through 3-6 the 


* 
fixed A technique provides data in an easily usable form, 


toa 








APPENDIX. B 


CHANGES TO COMPUTER PROGRAM IN REFERENCE 1 


The following changes were made to the computer program 
in Ref. 1 to convert to starred parameters. 
1. Program #1 
Statement number 0002 (COMPLEX *16 A,B) was deleted 
and two type declaration statements (REAL*8 TH) and 
(COMPLEX*16 A) were inserted. The namelist statement, 
number 0008, was revised to read: NAMELIST / LIST /N,REY, 
WReASTAR,PHIS;VEL. Statement number 0018, B = DCMPLX(BR,BI) 
was deleted and the following statements added: PHI = 
einsy 97.2958, AR = ASTAR * (DCOS(PHI)), AI = ASTAR *(DSIN 
Oo), THD.=, (TH*180.0)/3.141592654. Other changes to 
program #1 were those required to write out the revised 
Imomts . 
me ~‘oubToOuUtine. DETGEO 
One small change was made to this subroutine due to 
the fact that values required by external functions CHM1E1 
and CHM2E1l were passed by DEIGEO. Statement 0010 of DEIGEO, 
B = BETA, was deleted and TH = THETA was inserted. The 
common statement and type declaration statements were 
revised to incorporate the change to starred parameters. 
So. “bunetvons (Gn Wieland GHMZE 1 
External functions CHM1E1 and CHM2E1 required 


expensive MOGltication as follows: 
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oe A aT ne ee. 2 ee eo 


The type declaration statement (REAL*8 TH,DUR) was added. 
CH4M1(Y) = A/REY was changed to CH4M1(Y) = A*EI/REY. 
CH2M1(y) was changed to equal 
eee 2 PiDO-y**2)rZvO"ART*(A**2)/REY . 
CHOM1(Y) was changed to equal 
Seren (Ane 2 = (1. SPOeAET -CiDO- W4*2) - (A**2) /REY) 
Po0* At ) 
CH2M2(Y) = A changed to CH2M2(Y) = AEI. 
The following statements were added after CH2M2(Y) and 
ENTRY CHM2E1(k,Y): 
DUR = 0.0 
DU = DCMPLX(DUR, TH) 
EI @@GEXP (DU) 
AEI = Atel 


aeZ CDEXP(22DU) . 
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